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Amendments to the Claims: 

This listing of claims will replace all prior versions, and listings, of claims in the application: 
Listing of Claims: 

L-18. (Canceled). 

19. (Original) A method for detecting faults in a chiller based on vibration 
amplitude limits, comprising: 

calculating vibration amplitude limits of the chiller using statistics and 
historical data for the chiller; 

estimating an at least two-dimensional density estimate; and 

weighting the historical data based on when the historical data was generated; 

wherein the vibration amplitude limits are calculated as a function of 
frequency for an entire frequency spectrum. 

20. (Original) The method of Claim 19 further comprising removing outlier data. 

21. (Original) The method of Claim 20 wherein the at least two-dimensional 
density estimate utilizes frequency and amplitude directions of the frequency spectrum. 

22. (Original) The method of Claim 21 wherein the at least two-dimensional 
density estimate is a J-dimensional kernel density estimate. 
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23. (Original) The method of Claim 22 wherein the ^-dimensional kernel density 
estimate for point x of a dataset with n data points is given by: 

p(x) = ±-±\H\- m K(H-">(x-xj)) 

where, Xj is the 7 th observation of the dataset, K (w) is a ^-dimensional kernel, H is a 
bandwidth matrix, and | • | denotes a matrix determinant. 

24. (Original) The method of Claim 22 further including obtaining vibration 
spectra comprising individual spectrum for the chiller from a database. 

25. (Original) The method of Claim 24 further comprising calculating a frequency 
for the individual spectrum and identifying an individual spectrum having the smallest 
number of frequency lines. 

> 26. (Original) The method of Claim 25 further comprising calculating noise 

bandwidths and a largest noise bandwidth. 

27. (Original) The method of Claim 26 further comprising collecting vibration 
data from all spectra in a given frequency range. 

28. (Original) The method of Claim 19 further comprising calculating a 
conditional kernel density. 

29. (Original) The method of Claim 28 wherein calculating the conditional kernel 
density comprises estimating an unknown probability density for a given dataset. 
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30. (Original) A method for determining vibration amplitude limits of a 
mechanical device comprising: 



identifying a mechanical device and a frequency range for a spectrum to be 



retrieving vibration spectra comprising individual spectrum for the mechanical 



calculating frequency for the individual spectrum; 

identifying the individual spectrum with a smallest number of frequency lines; 

calculating noise bandwidths and a largest noise bandwidth; 

removing outlier data; 

calculating conditional kernel density; and 

calculating vibration amplitude limits to detect faults in the mechanical device. 



31. (Original) The method of Claim 30 wherein the mechanical device comprises 
a chiller for an HVAC system. 

32. (Original) The method of Claim 30 wherein the vibration spectra for the 



33. (Original) The method of Claim 32 wherein calculating conditional kernel 
density comprises estimating an unknown probability density for a given dataset. 

34. (Original) The method of Claim 33 wherein the probability density estimate at 
a point x for a one-dimensional dataset with n data points is given by: 



where, xj is the j observation of the dataset, h is a bandwidth that characterizes a spread of 
the kernel, and k (•) is a kernel density function that is symmetric and satisfies the condition: 



analyzed; 



device and the frequency range; 



mechanical device and the frequency range is obtained from a database. 
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35. (Original) The method of Claim 33 wherein the kernel density estimate is at 
least a two-dimensional kernel density estimate utilizing frequency and amplitude directions 
of the frequency spectrum. 

36. (Original) The method of Claim 35 wherein a rf-dimensional kernel density 
estimate is given by: 

p(x) = l-£\H\- U2 K(H-v\x-xj)) 

where K (w) is a rf-dimensional kernel, H is a bandwidth matrix, and | • | denotes a matrix 
determinant. 

37. (New) A method for determining vibration amplitude limits to detect faults in 
mechanical equipment, comprising: 

estimating a data probability distribution based on data for the mechanical 

equipment; 

utilizing the data probability distribution to calculate the vibration amplitude 
limits, such that the data probability distribution is calculated using statistics and historical 
data of the mechanical equipment including using a kernel density method, wherein the kernel 
density method comprises calculating conditional kernel density including estimating an 
unknown probability density for a given dataset, the probability density estimate at a point x 
for a one-dimensional dataset with n data points being given by: 

where, xj is the / h observation of dataset X, h is a bandwidth that characterizes a spread of the 
kernel, and k (•) is a kernel density function that is symmetric and satisfies the condition: 

^K{u)du =1; 

removing outlier data; and 

calculating the vibration amplitude limits as a function of frequency for a 
substantial portion of the frequency spectrum. 
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38. (New) The method of Claim 37 wherein the kernel density estimate is a two- 
dimensional kernel density estimate utilizing frequency and amplitude directions of the 
frequency spectrum. 

39. (New) The method of Claim 38 wherein a ^-dimensional kernel density 
estimate is generally written as: 

p{x) = ±±\H\ V2 YL(H-»\x-xA) 

where K (u) is a c/-dimensional kernel, H is a bandwidth matrix, and | ■ | denotes a matrix 
determinant. 
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